Controlling Strong Electromagnetic Fields at a Sub- Wavelength Scale 
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We investigate the optical response of two sub- wavelength grooves on a metallic screen, separated 
by a sub-wavelength distance. We show that the Fabry-Perot-like mode, already observed in one- 
dimensional periodic gratings and known for a single slit, splits into two resonances in our system : a 
symmetrical mode with a small Q-factor, and an antisymmetric one which leads to a much stronger 
light enhancement. This behavior results from the near-field coupling of the grooves. Moreover, the 
use of a second incident wave allows to control the localization of the photons in the groove of our 
choice, depending on the phase difference between the two incident waves. The system exactly acts 
as a sub-wavelength optical switch operated from far-field. 

PACS numbers: 71. 36-l-c,73.20.Mf, 78.66. Bz 
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Surface Enhancement Raman Scattering (SERS) still 
remains a mystery in a large part, even though it is now 
accepted that the excitation of localized electromagnetic 
modes of irregular metallic surfaces is involved in its ba- 
sic mechanism0, Optical excitation of such modes 
can indeed lead to important concentration of electro- 
magnetic energy in volumes (cavities) much smaller than 
A'^ where A is the excitation wavelength, as it is the case 
for SERS active surfaces. These specific places of very 
strong electromagnetic fields localization are called "ac- 
tive sites" or "hot spots". However, the debate on the 
origin of these hot spots remains open, as well as the 
hope to control one day this phenomenon. The large in- 
terest raised by this fundamental physics is also increased 
by its wide potential applications in biochips, sensors, 
nano-antennae, optoelectronics or energy transport on 
nanostructured surfaces. 

In this letter, we consider a simple system which al- 
lows to produce and control the localization in space of 
such hot spot phenomenon. It only consists of two deep 
grooves on a plane metallic gold surface (fig.l). The ex- 
cited modes appear, for the chosen geometry, in the in- 
frared region where we can consider the metal as be- 
ing a good reflector. Under this condition, a reliable 
theoretical method, i.e the modal method using surface 
impedance boundary conditions, can be used [3j. This 
method has already demonstrated its ability to give a 
good qualitative and quantitative agreement with the 
measured reflectivity of metallic gratings 0, |H 01 ■ The 
case of one groove only was considered a long time ago 
0, while the transmission for one 8j and two slits Q 
were only recently considered. In contrast with 0, the 
distance between our two grooves is small with respect to 
the incident wavelength. Very recently it also was shown 
|lO| that sharp and deep resonances appear in the trans- 
mission response of gratings with more than one slit per 
period or in gold dipole antennas [fl). We here analyze 
the physical origin of this new kind of resonances for a 



two slit system. As we will see, this allows us to point 
out some very fundamental aspects of electromagnetic 
resonances on metallic surfaces, and to control the light 
localization by using a simple device. 

We consider a p-polarized incident plane wave (elec- 
tric field in the plane of incidence) with a wavevector 
k — 2tt/\ impinging on the surface at an angle 9 (fig.l). 
The knowledge of the magnetic field in the z-direction 
completely solves the problem as — Hy = = 0, 
Ex = {i/ckeo)dHz/dy and Ey — {—i/ckeo)dHz/dx. In 
region (I) , the field is expressed as the sum of the incident 
wave and the reflected ones by: 
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where the distribution R{Q) represents the amplitude 
of the reflected field at the wavevector {Q, q) with q = 
■\/fc2 — Q^. In region (II) one has: 
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where Ii{x) (respectively l2{x)) equal 1 in the interval 
[{-w-d)/2,{w-d)/2\ (resp. [{d-w)/2,{d + w)/2]) imd 
zero elsewhere, a = [(1 Z)/(l - Z)]$^ with $ = e"*'''', 
Z — l/\/e is the surface impedance of the metal and e 
its dielectric constant. The expression for iji^^^ assumes 
that the field is constant along x within each groove, 
which is a good approximation in the limit where w << X 
0, 01 . To illustrate our results numerically, we have fixed 
w = 0.2 fj,m, h = 1.5 /im, and d = 0.5 fim all along the 
paper. The values of the complex dielectric constant e(A) 
are taken from ^3 • 

The unknown variables are the distribution R{Q) and 
the field amplitudes Ai and A2 in the first and sec- 
ond groove respectively. A set of equations is obtained 
by applying the boundary conditions at the interface 



2 










k .® H^.e''* 


-9 










FIG. 1: Geometrical configuration and parameters. Region 
(I) and (II) respectively corresponds to the region above and 
below the metallic surface. is the phase difference between 
the two incident waves described at the end of the paper. 



2/ = : hS!'^ — H^^^^ at the mouth of each groove, 
and dnP/dy + ikZHi'^ = dH^"^ /dy + ikZHi^^'> along 
the whole interface. After some elementary algebra (see 
for detailed procedm-e), the vector A = (^1,^2) 
is related to the excitation vector V = (14, V2) (null 
without the incoming wave), by the matricial relation 
A = M~^V, where M is the 2x2 symmetrical matrix 
which verifies mn — 77122 with: 
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-00 q 



dQ, 



q + kZ 

where F = w/X. The coordinates of the vector V are: 



2 cos 6* 
cos 6* + - 



and V2 = e"^Vo with: 
- sec (fc sm{9)w/2) 



where we have introduced the angle ip = kdsin{6)/2. 
The matrix M has two eigenvalues ■m± = — $^)e± 
with respective eigenvectors XJ± = (1,±1) and: 

e± = ■:^^-^{cot{kh)-iZ)~2ir{l + Z)x 
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We have made the variable change Q — ku in the inte- 
grals. The solution of the problem is then: 
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FIG. 2: _E_Fi of the cavity centered at a; = —d/2 calculated 
for = 0° (full line) and 6 — 30° (dotted line) as a function 
of the wavenumber. Inset a) represents the behavior of EF2 
(cavity centered at 2; = +d/2), as a function of the incidence 
angle 6 and at 1484 cm^^. Inset b) gives EFi and EF2 of 
both cavities calculated for 50°. 



At the sight of eq.(l), one can see that the system 
presents two electromagnetic resonances at fc = k±, 
which appear when 3fi(e+) = and 5R(e_) = 0, with line- 
shapes respectively governed by 3(64.) and 3(e_) (3fi(x) 
and 9(a;) being the real and imaginary parts of x). The 
fields in the cavities are always a linear combination of 
the two eigenvectors A ~ a_U_ -I- a-|_U+. However, 
when k — k+ (respectively k = k^), the vector A is al- 
most coUinear with U_|_ (resp. U_) as the amplitudes in 
the two cavities are dominated by the same (resp. the 
opposite) term. We will thus call the resonance occur- 
ring at k — k^ the (— ) antisymmetric mode and that 
occurring at /c = fc+ the (-I-) symmetrical one. Contrary 
to the (-I-) mode which always exists, the (— ) one only 
develops for 6^0 (see fig. 2) as it vanishes at normal 
incidence with tamp = 0. Its bandwidth is much thinner 
than that of the symmetrical mode and its enhancement 
factor is much larger. The enhancement factor {EF), de- 
fined as \Ex/Eq\'^ where Eq is the incident electric field, 
refiects the amount of stocked electromagnetic energy at 
the resonances. For convenience, we note EFi and EF2 
the enhancement factors calculated at the mouth of each 
cavity, i.e at a; = ±c?/2 and y = where they are ex- 
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pressed as EFn=i.2 = |^r«=i,2(l - The EF of the 
(— ) mode, shown in the inset (a) of fig. 2, increases with 
6 and its value can reach more than 10^ whereas that of 
the symmetrical mode stays at around 100. 
Another important point to highlight is that around the 
(— ) resonance, the fields in the two cavities are not 
strictly identical. Inset (b) of fig. 2 displays the EF at 
the mouth of each cavity close to A; = At 1483 cto~^ 
the maximum of EFi is reached whereas the value of 
EF2 is still low; at 1490 cm'^ EFi and EF2 both take 
the same value. Around this mode, the system thus de- 
velops a very high sensitivity: with a very little variation 
of wavenumber (here less than one percent), the field 
"jumps" from one cavity to the other. This behavior is 
qualitatively comparable to the "unstable" behavior of 
hot spots observed on SERS active surfaces. 

In the following, we consider the metal as being a per- 
fect reflector, i.e Z = 0. This approximation induces only 
small quantitative modifications and allows an analytical 
study which highly helps to clarify the physics of the 
problem, however, the presented numerical results are 
obtained without this approximation, i.e using the finite 
value of e(A). We first compare the two grooves system to 
the one where there is only one groove centered at x = 0. 
In this case, the amplitude of the field Aq in the unique 
cavity is given by ^0 = * [l — ^'^] K)/e, with : 



e = cot{kh) - 2iT 



+00 



sec^{kwu/2) 



du. 



The resonance of this cavity occurs at k — — luq/c 
for which 5R(e) = 0. Close to fcp, the field ^0 can be 
expanded around loq as: 



^0 



Co 



luq-lu - «7o/2' 



with C'o ~ icVo/2h, and where we have taken advan- 
tage of the fact that at the resonance koh w 7r/2 [^. 
This equation is typical of a forced oscillator and, as the 
electric field inside the cavity is proportional to Aq , indi- 
cates that the cavity behaves as a forced oscillating dipole 
with a radiation damping 70 = 2wluq /ttc and an effective 
electromagnetic radius tq — 2w/tt. The effective dipo- 
lar momentum, parallel to the interface, takes its maxi- 
mum at the mouth of the groove and decreases along the 
vertical walls. The maximum of intensity at uj — luq is 
l^oP ~ 4/(/cow)^, typically of order 100 for our geomet- 
rical parameters. We now expand, in the same manner, 
the values of e± around the same kg for the two groove 
system. We easily get: 



— lu — /2) h/c 
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FIG. 3: EFi calculated for 9 — 20° showing the resonances 
and (— ) characterized by their eigenfrequencies uj±, lo- 
cated on both sides of the frequency resonance uiq of a unique 
cavity, and their bandwidth y± (a). Right part schematically 
represents the in-phase coupling of the mode (b) and the 
anti-phase coupling of the ( — ) mode (c) and their correspond- 
ing equivalent dipole. 



shift A, of the order of 70 <^ wo, is: 

^ 7o Z"^"" cos (fcodw) sec^ (fcowu/2) 



1 



du 



(3) 



with ll!± = Wo T a, 7_|_ = 270 and 7_ = 70 {kod/2) . The 



Eq. (2, 3) confirm our numerical observation as they 
show that the width of the (— ) mode, driven by 7_, is 
much lower than that of the mode, driven by 7+, 
owing to the small factor (kod)"^ (and recalling our sub- 
wavelength coupling hypothesis : Ao >> rf). A phys- 
ical image of these resonances can be given noticing 
that our results are completely similar to those obtained 
by Lvuboshitz|l3l| for two near- field coupled oscillating 
dipoles. Our resonances thus arise from the near-field 
coupling of two identical grooves, individually resonat- 
ing at ujQ. The symmetrical (-I-) mode corresponds to 
the in-phase oscillation of each cavity whereas the sec- 
ond one corresponds to an anti-phase oscillation. The 
distribution of electric field in the cavities for each mode 
is sketched in fig. 3. As a consequence of this coupling, 
the (-I-) mode has a strong dipolar character with an ef- 
fective dipolar moment close to twice that of a unique 
cavity and a large electromagnetic radius r+ = 2ro. On 
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the opposite, the (— ) mode has an effective dipolar mo- 
ment almost null, with a much smaller electromagnetic 
radius r_ = ro{kod/2)^, and its radiation pattern is es- 
sentially that of a quadripole. This explains why this 
mode is weakly radiative and with an extremely narrow 
lineshape, very different from the width of the in-phase 
mode. 

Searching for the location of the maximum of the field in 
each cavity around the (— ) mode, ones gets for non zero 
6: 



,n=l,2 



U!- — 



(-ir(¥)So 



4sin6' I 1 
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+ O{{kodf) 



{kowy{kod) 
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where l^maxP is proportional to the intensity of the field 



in both cavities at lu„ 



The two maxima 



and 



^max separated by a very small frequency difference 
of the order of (fcorf)^7o, which, together with the narrow 
lineshape of the resonance, explains why the profile of 
the field strongly varies in this region. The magnitude 
of |j4TOaxP requires some comment. Indeed, for an usual 
oscillator with damping 7, the maximum of intensity of 
the oscillation scales as 7"-^, so that should scale 

as 7!^ ^ {kod)~^ instead of (fcod)~^. The field intensity 
of the (+) mode scales, as expected, as 7^^ (eq. 1). The 
reason for that is that the (+) and (— ) modes are not 
sensitive to the same parts of the incident electric field. 
Since d/X <C 1, the latter can express at interface as 
£"0(1 + ikx) at the scale of our two-grooves system. The 
even term corresponding to the mean value of the field 
excites the (+) mode and the odd one, corresponding to 
the local variations of the field, excites the (— ) mode. 
This mode is thus sensitive to an "effective" field of in- 
tensity ~ £'□(^0^)^ at the mouth of the grooves, whereas 
the mode is excited by an effective field of intensity 
Eq. This is the origin of the lost of a factor (fcorf)"^ in 
the intensities of the mode (— ). The latter results from 
a strong resonator, but excited by a very weak effective 
field. 

We now take advantage of our understanding to control 
- from far field - the light localization in the cavity of our 
choice, or in both. To do so, we introduce a new free pa- 
rameter by sending a second incident plane wave, at the 
same frequency, with an incidence angle —9, and with a 
phase difference (p with respect to the first incident wave 
(fig.l). Changing (f), we can control the incident effective 
fields respectively exciting one mode or the other. Dif- 
ferent states, that we code as : (1,1), (1,-1), (1,0) and 
(0, —1) can be achieved. The first two, (1, 1) and (1, —1) 
respectively correspond to the case where only the pure 
{+) or only the pure (— ) resonances are excited. The 
cavities arc then completely in-phase or in anti-phase. 
The two other ones correspond to cases where one of the 
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FIG. 4: EFi and EF2 spectra at 9 — 45° for two incoming 
waves (a) and related mappings of the electric field amplitude 
Ex (b, c) at 1490 cm^^. For = 0(1.0) tlie first cavity is lit, 
EFi is the dotted line, and the second cavity is extinguished, 
EF2 is the full line. For <p = (p(o,-i) it is the opposite. 



cavities is lit (cavity 1 for (1, 0), and cavity 2 for (0, —1)). 
As (/) is a parameter easy to modify, for instance chang- 
ing the optical path, we can control in straightforward 
manner the field localization. 
With two incoming waves, the field becomes: 



rr. Aksin$x i A 

J^in.r.. — e ~j~ e 



-ksinOx) ^—ik COS Oy 
e- , 



and the solution for each cavity can be written as: 

cos{4>/2) . sin((/)/2) 
e+ e_ 

where we did not write explicitly some unimportant pref- 

actor common to both cavities. From these equations, 
it is easy to see that for </> = (t>(i^i) = 0, one gets 
Ai= A2 ^ 1/e-i-, so that at fc = fc+ we have the pure (-I-) 
resonance. In the same manner, the pure (— ) resonance 



can be excited at k = k- when 



-1) 



TT where 



Ai = —A2 ~ — tan(^/e_. This last state presents an 

extremely high EF - 10'' at 6* = 80°. 

More subtle is the possibility to control the extinction 

of the field in only one of the cavities (of our choice) 

while the other one is resonating. Eq.(4) shows that this 
can be achieved provided that cot(/)/2 = ±(e+/e_) tantp. 
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the sign " + " corresponding to the (0,-1) state and the 
" — " sign to the (1,0) state. This condition can be sat- 
isfied (since the function cot(x) can vary from — cxd to 
+oo), provided that e+/e_ is real. This is obtained for 
« + 2A(7_/7+), which is very close to a;_. Fig- 
ure 4 represents the EF of both cavities either choosing 
(j) = 4>(ifi) or = 0(o^_i), together with the related map- 
pings of the electric field amplitude E^- These show how 
the field can be strongly localized in only one of the cav- 
ity, while the second one is completely extinguished and 
this, even though the cavities are identical and separated 
by a sub-wavelength distance. 

In conclusion, we have demonstrated that the near-field 
coupling of two metallic resonating cavities leads to a res- 
onance with an extremely thin spectral width, which can 
localize very intense fields. This could be a key point in 
the understanding of the SERS, as the described physics 
should remain valid in the visible region, except for a 
scaling factor. Finally, we proposed a simple way to con- 
trol the near-field of each cavity, enabling this system to 
act as a sub-wavelength optical switch simply operated 
from the far-field. 

This work is partly the result of illuminating discus- 
sions of P. Quemerais with D. Mayou which have been 
very beneficial to us. We also would like to thank P. 



Lalanne for helpful conversations about electromagnetic 
resonances in gratings. 
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